TUM-HEP 871/12 
FLAVOUR(267104)-ERC-33 
KIAS-PREPRINT-P 12084 
SLAC-PUB-15323 



B — y D(*>tv and B — >• rv in chiral U(l)' models 
with flavored multi Higgs doublets 

P. Ko, 1 Yuji Omura, 2 and Chaehyun Yu 1 ' 3 

1 School of Physics, KIAS, Seoul 130-722, Korea 
2 Physik Department T30, Technische Universitat Miinchen, 
James-Franck-Straj3 e, 85748 Garching, Germany 

3 SLAC National Accelerator Laboratory, 
2575 Sand Hill Rd, Menlo Park, CA 94025, USA 
(Dated: December 20, 2012) 

Abstract 

We discuss semileptonic and leptonic B decays, B — > D^*'tu and B — > tu, in the chiral U(l)' 
models which were proposed by the present authors in the context of the top forward-backward 
asymmetry (^4 FB ) observed at the Tevatron. In these models, extra Higgs doublets with nonzero 
J7(l)' charges are required in order to make the realistic mass matrix for up-type quarks. Then the 
extra (pseudo)scalars contribute to ^4 FB with large flavor-changing Yukawa couplings involving top 
quark. The contribution of the charged Higgs to Ap B is negligible, but it may significantly affect 
B decays: especially, B — > D^*'tv and B — > tu. We investigate constraints on the B decays, based 
on the recent results in BaBar and Belle experiments, and discuss the possibility that the allowed 
parameter region in the B decays can achieve large ^4p B - 
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I. INTRODUCTION 



Recently, the BaBar collaboration analyzed semileptonic B decays, B — > Dtu and B — > 
D*tu [1], and investigated the ratios of the branching ratio for B — > D^*'rv to B — > D^*Hu 
(I = e,/x), 

#(£>(*)) = BR(B D [ * ] tv) /BR(£ D^lu). (1) 

The results are R(D) = 0.440 ± 0.072 and R(D*) = 0.332 ± 0.030 which deviate from the 
Standard Model (SM) predictions, R(D) SM = 0.297± 0.017 [2] and R(D*) SM = 0.252 ±0.003 
[3] by 2.2a and 2.7cr, respectively [4]. The combined discrepancy is about 3.4a [1], which 
might be an evidence of new physics as discussed in some recent works [4-10]. One good 
candidate for new physics for this anomaly is a charged Higgs boson in the extended SM 
with extra Higgs doublets [4, 5]. 

On the other hand, a leptonic B decay B — > tv was measured at BaBar [11] and Belle [12], 
the average of which is BR(5 tv) = (1.67 ± 0.3) x 10~ 4 [13]. In the SM calculation, 
there are some uncertainties from \V u b\ and the B meson decay constant Jb- Still the 
measured number is slightly inconsistent with the SM prediction, for example, given by 
UTfit Collaboration, BR(B rz/) SM = (0.84 ±0.11) x 10~ 4 [14]. The Belle experiment 
recently presented a new result, BR(5 — > tv) = 0.72^Q25±0.11 by making use of a hadronic 
tagging method for r decays with the full data sample [15], and the combined average is 
consistent with the SM prediction. Since both semileptonic and leptonic B decays may be 
affected by the same new physics (e.g., charged Higgs boson), such new physics scenario will 
be strongly constrained by combined analysis of the B decays. 

An interesting point is that it is difficult for the so-called type-II two-Higgs-doublet 
model (2HDM), which could be well motivated by minimal supersymmetric SM (MSSM), to 
explain the discrepancies in R(D) and R(D*) [1]. Other types of 2HDMs with natural flavor 
conservation, where Yukawa couplings are controlled by Z 2 [16] or U(1)h [17] symmetry, 
also allow only the so-called minimal flavor violation (MFV) in the charged Higgs sector, 
and it is impossible to accommodate R{D) and R(D*) at BaBar simultaneously [4, 5, 9, 10]. 
Eventually, we would have to introduce non-minimal-fiavor-violating (non-MFV) terms in 
order to explain the R(D) and R(D*) data, which would tend to generate too large flavor 
changing neutral currents (FCNCs). It would be highly nontrivial to introduce non-MFV 
interactions in the Yukawa couplings in 2HDMs without too excessive flavor violations in 

2 



the K and B meson sectors. 

In Refs. [18], the present authors proposed flavor models with extra Higgs fields, where 
gauged U(l) ; symmetry forbids the potentially problematic FCNCs for B-B and K-K mix- 
ings but allows certain amounts of FCNCs which are still consistent with experimental 
data, by slightly breaking the criteria of Ref. [17]. There, neutral CP-even scalars and a 
CP-odd scalar can have large (t, q) q=u c t elements of Yukawa couplings. Their tree-level me- 
diations enhance top forward-backward asymmetry (Ap B ) at the Tevatron, accommodating 
with the newest strong constraints from LHC, thanks to the destructive interference among 
the scalars [18, 19]. In previous works, phenomenology of the charged Higgs boson in our 
models was not considered carefully because the charge Higgs boson cannot have sizable 
contribution to the top quark production at hadron colliders. However, couplings of the 
charged Higgs boson to the bottom quark may be large so that the models could be strongly 
constrained by the B decays. For example, the (b, u) element of the charged Higgs coupling, 
which is constrained by B — > tu, may become large if the (t, u) element of the Yukawa 
coupling of the pseudoscalar boson is large. Besides, the (b, c) element, which modifies the 
branching ratios of B — >■ D^'ru, could be large, if the (t, c) element of the pseudoscalar 
Yukawa coupling is large (see Eq.s (8), for example). In this paper, we first investigate if 
our U(l)' flavor models can explain the discrepancies in the B semileptonic decays while 
keeping consistency with BK(B — > tv). Then we discuss the possibility that the allowed 
parameter regions can achieve the large enough Ap B which were observed at the Tevatron. 

This paper is organized as follows. In Sec. II, we give a brief review of our models proposed 
in Ref. [18]. We assign family dependent U(l)' charges to the right-handed (RH) up-type 
quarks in order to generate flavor changing Z' — ur — couplings in the mass eigenstates. 
Then we will immediately realize that it is mandatory to introduce extra U(l)'- charged 
Higgs doublets in order to write Yukawa couplings for the up-type quarks including top 
quark, which have been first realized in Refs. [18]. Thus we are led to multi-Higgs doublet 
models in the presence of a new spin-1 Z' boson with chiral couplings to the SM fermions. 
Sections III and IV are devoted to discussion of the phenomenology for B — > D^*<tv and 
B — > tv in our 2HDM and three-Higgs-doublet model (3HDM), respectively. We also give 
comments on the constraints on the charged Higgs boson from the B — > AT, 7 process in 
Sec. V. Finally, we summarize our results in Sec. VI. The general Higgs potential in multi 
Higgs doublet models is discussed in Appendix A. 
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II. MODELS WITH EXTRA HIGGS AND GAUGED FLAVOR 17(1)' 



Adding extra Higgs doublets to the SM is one of the interesting extensions to the SM. 
However such an extension generically suffer from neutral Higgs-mediated FCNC problem at 
tree level, if up- and down-type quark masses get contributions from the vacuum expectation 
values (VEVs) of more than one Higgs doublets. For example, when the mass matrix of up- 
type quarks (M£) depends on two Higgs doublets, (Hi,H 2 ), as 

(M u ) lj = (y 1 ) i ,(H 1 ) + (y 2 ) lJ (H 2 ), 

all scalar components of (Hi,H 2 ) would have flavor-dependent couplings. This is because 
{yi)ij and {y 2 )ij cannot be diagonalized simultaneously without any relation between (y\)ij 
and (y 2 )ij- A simple way to control the flavor structures of Yukawa couplings is to assign 
a symmetry to extra Higgs and matter fields. The most popular symmetry is an extra 
Z 2 symmetry which is softly broken [16]. In Refs. [17] and [18], gauged U(l)' symmetry 
is assigned to the SM fermions and newly introduced extra Higgs doublets instead of Z 2 
symmetry. Especially, in Refs. [18], only right-handed (RH) up-type quarks are charged 
flavor-dependently under U(l)', and FCNCs involving only top quark can be enhanced. The 
authors constructed both 2HDM and 3HDM depending on the U(l)' charge assignments to 
the SM fermions. 

When only the RH up-type quarks are charged in a flavor dependent way, the extra Higgs 
doublets charged under U{1)' couple with the RH up-type quarks in the form 

V y = ytQJIjUjv + vtjQi^DRj + y\^H 2 E Rj + h.c. (2) 

Here Hj(j = u, c, t) are charged under U (1)', but H 2 is neutral under U (1)' and has the same 
quantum numbers as the SM Higgs doublet. The SM leptons and the down- type quarks get 
masses from VEV of H 2 . Note that the up- type quarks cannot have gauge invariant mass 
terms without new £7(l)'-charged Higgs doublets Hj =UtCtt . 

In general, there may be up to four Higgs doublets: H 2 and H uct . The actual number 
of Higgs doublets will depend on the U{\)' charge assignment. Motivated by the Ap B at 
the Tevatron, we use the charge assignment (uj) = (0,0, 1) on the RH up-type quarks in 
the 2HDM, where we identify H u and H c as H 2 and H t as Hi, respectively. In the 3HDM, 
we use the charge assignment (uj) = (—1,0, 1) and we identify H u as Hi, H c as H 2 , and 
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H t as if 3, respectively *. In addition to the extra Higgs doublets Hj, a SM singlet scalar 
with nonzero U(l)' charge, $, may exist in order to break U(l)' and mix with CP-even 
components of Hj, after electroweak and U(l)' symmetry breaking. In the following two 
sections, we will concentrate on these 2HDM and 3HDM. Then we investigate the impact 
of the charged Higgs boson sector on the (semi)leptonic B decays and discuss the results in 
the context of the A^-n at the Tevatron. 



III. 2HDM 



In this section, we consider the chiral U(l)' flavor model with the charge assignment 
= (0,0,1) in the interaction eigenstatesL There are a neutral Higgs doublet, H 2 and 
U(l) '-charged Higgs doublet, Hi, as we discussed in Sec. II. 



A. Yukawa couplings 

In the 2HDM case, there are three CP-even neutral scalars (hi, h 2 , h 3 ), one CP-odd 
pseudoscalar (a), and one charged Higgs pair (h^) after electroweak and U(l)' symmetry 
breaking. Note that there is one more CP-even scalar from the £7(l)'-charged singlet scalar 
$ compared with the usual 2HDM. 

Let us define the Yukawa couplings for the neutral scalar bosons in the mass bases as 
follows: 

v v u = Y ij ik) hkuTiURj - iY^auTiURj + h.c. (3) 
yu(k),au ( j e p en( j on Higgs VEVs and the diagonalizing matrices for the quark mass 
matrices. The Yukawa couplings can be derived in a straightforward manner. For example, 
the Yukawa couplings of the lightest CP-even scalar boson hi and pseudo scalar boson a are 
given by 

u(l) Tfl^ COS (X ^Tfl^ 

Y ii = — 7T cosa^ij H _ L-^^ . s i n ( a _ ft) cosa$, (4) 

J vcosp vsmzp 



m" tan (3 2m" 

Ojj 

v v sin 2/3 

The parameters v and (3 are defined by ((Hi), {H^}) = (vcos/3/\/ r 2,vsm/3/\/ r 2), and a and 

a$ are the mixing parameters among the 3 CP-even neutral scalars. 

* There could be other assignment of (v,j) but we consider only these two cases for simplicity, 
t We use the same notations as in Refs. [18]. 



In the above equations, the mixing matrix g R is defined as [18] 

{glh = {qDI = (RuhkU k (R u )* k , (6) 

where the matrix (R u )ij is defined by 

(AfjMOo- = (R u )i(m 2 u ) k (R u ) kj . 
For the 2HDM case with (u k ) = (0, 0, 1), g R is reduced to 

{9b)%3 = (dnTji = (Ru)i3(Ru)* j3 - 

Note that the neutral Higgs bosons have flavor-dependent couplings to the up-type quarks 
(see the second terms in Eqs. (4) and (5)). The flavor-dependent couplings of Z' which is 
the gauge boson of U(l)' are also linear in {g R )ij in our 2HDM [18]. These flavor-dependent 
couplings make additional contributions to Ap B at the Tevatron, tt cross section, and the 
same-sign top-quark pair production cross section. Recently, the CMS collaboration has 
announced stringent bounds on the cross section for the same-sign top-quark pair production: 
a tl < 0.39 pb at 95% confidence level [20, 21]. This strong bound excludes simple scenarios 
such as the original Z 1 model that only one mediator contribution is taken into account. 
However, flavor models usually have several mediators to couple with the SM particles 
flavor-dependently, and they interfere with each other. In fact, our pseudoscalar, CP-even 
scalars, and Z' have destructive interference in our U(l)' flavor models and could find the 
points evading the stringent upper bound on the same-sign top-quark pair production and 
enhancing Ap B [18, 19]. 

The origin of this non-MFV couplings in our models is the flavor-dependent U(l)' inter- 
actions of Z'. (Recall that (g^)^ oc Sij if the U(l)' were flavor-independent (u k ) oc (1, 1, 1).) 
Thus our 2HDM is a nice realization of non-MFV models where the flavor non-universality 
has its origin in the flavor-dependent U(l)' gauge interactions. The amount of flavor non- 
universality is related with the local gauge symmetry and its spontaneous breaking, and 
is not completely arbitrary or ad hoc. These are very unique features of the 2HDM (and 
3HDM described in the next section) we have proposed in Refs. [18]. 

Similarly to the neutral Higgs Yukawa couplings, the charged Higgs Yukawa couplings 
are defined by 

V c = -Y^ u hT~d^u R j + Y^ d h + uUd Rj + Y^ l h + VTie Rj + h.c. (7) 
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There are definite relations between Yukawa couplings of the pseudoscalar boson (a) and 
those of charged Higgs bosons (h 1 ^): 

i 

3^ = £^ y «V2> (8) 
i 

I 

where V is the CKM matrix. 

In order to accommodate the large Ap B observed at the Tevatron, large flavor-changing 
Yukawa couplings of the pseudoscalar are inevitable [18]. As the authors pointed out in 
Refs. [18, 19], the pseudoscalar contribution to the same-sign top-quark pair production has 
a destructive interference with other neutral scalar and Z' contributions. Once we consider 
a large FCNC in the pseudoscalar Yukawa couplings in the up-type quark sector in order 
to enhance the Ap B , we will necessarily have large flavor-changing couplings in the charged 
Higgs sector, which would affect various phenomenology including the B meson system. 

In the following, we will discuss phenomenology of the charged Higgs boson. 

B. «(£>(*>) and BR(B ->■ tv) in 2HDM 

The charged Higgs boson in our model will contribute to the (semi)leptonic B decays 
and would modify the SM lepton universality which is a result of the W contributions 
derived from the underlying SU {2) L gauge theory. Since the charged Higgs contributions are 
proportional to the final lepton mass, we consider the contributions only to B — y (D^)ru 
and B — y tv, and use the SM predictions for other leptonic channels. 

In the effective Hamiltonian approach, the semileptonic decays B — y qru and leptonic 
decay B —y tv are described by the effective Hamiltonian [4], 

#cff = Cf M {Wl,b L ){m^L) + Ct{qEb R ){r^v L ) + ct(<mb L )(w L ). (9) 

where q = u, c is the up-type quark flavor. In the above equation, is the Wilson 
coefficient for the W exchange in the SM and Cf! L are those for the charged Higgs exchange 
present in our mo dels. R(D^) and BR(B tv) are given by the following expressions 
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c c R b + cf 



depending on the above Wilson coefficients [4]: 
R(D) = R SM (l + 1-5 Re 

R(D*) = 



SM 



rich rich 

1 + 0.12 Re' R L 



+ 



C j* + Cf 



ricb 



rich 



0.05 



fjcb rjcb 

R L 



rich 



BR(B -> tv) 



r 2 it/ i 2 



57T 



m T f B m B T B 1 ^ 



?72 



2 \ 2 



B 



mbm 7 



Cr-C 

riub 



ub 



(10) 

(11) 

(12) 



where each Wilson coefficient is at the B meson scale [6]. Here, R^m are given by BK(B — > 
D^tv)/BR(B -> D^lv) in the SM. 

Integrating out the heavy degrees of freedoms (h^) in our 2HDM, we can obtain the 
Wilson coefficients at the charged Higgs mass scale: 



c 



qb 



ftqb 



m ^tan 2 /3 



Vib rrtfm T {g\ 



m 



h+ 
rribm T 
ml. 



RJk 

^ V qb m 2 h+ cos 2 ' 



tan 2 p. 



(13) 
(14) 

(15) 



We note that C R /C^ M is flavor-blind, but C^/Cq M depends on the flavor, q = u or c. If 
(dfdij = $ij were satisfied, Cjf R would be the same as the type-II 2HDM and one could not 
accommodate the R{D) and R(D*) simultaneously. However in our flavor- dependent U(l)' 
models, we have (g R )ij ^ $ij an d there appears a new term in Eq. (14) which is absent in 
the Type-II 2HDM. This new term may give rise to a possibility to accommodate R{D) 
and R{D*) unlike the Type-II 2HDMs. As discussed in the previous subsection, we could 
generate this non-MFV interactions from flavor-dependent U(l)' gauge couplings, which is a 
very interesting aspect of our multi-Higgs doublet models with flavor- dependent U(l)' gauge 
interactions. 

The contribution of the charged Higgs boson to the semileptonic B decays, B — > D^tv, 
is controlled by C L b R while that to the leptonic B decay, B — > tv, is affected by C^ b R . As 
discussed in Refs. [1, 4-7], it is impossible to explain R(D^) and BK(B — > tv) simultane- 
ously within the 2HDMs with MFV, where the Yukawa couplings are fixed by the angles 
a and ft up to the quark mass. However, if C c L h R and C L b R are mutually independent as in 
our models, both the measured values of R(D^) and BK(B — > tv) might be explained by 
taking appropriate parameters for the Wilson coefficients [4]. 
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In Refs. [18], the present authors suggested that large is required in order to achieve 
the large A^ B and evade the strong bound from the same-sign top-quark pair production 
signal. Such large can be realized by large (g R )tu and small sin 2/3, but it leads to a 
large Y fe ~ u according to Eq. (8) in our 2HDM. Therefore the charged Higgs contributions to 
the (semi)leptonic B decays could be too large, and have to be studied carefully. 

In numerical analysis for (semi)leptonic B decays in our 2HDMs, we take the following 
parameter regions: 1 < tan/3 < 100, 200 GeV < m h + < 1 TeV, and < \(g u R )tu\, \(g u R )tc\ < 1, 
respectively. We use the input parameters in the SM which are given by the global fit [22]. 

First, we consider only the B — > rv decay. In Fig. 1 (a), we show constraints on 
\C£ /CsmI an d I^rV^smI- The re d points are allowed by the combined data by Heavy 
Flavor Average Group (HFAG) [13] without the recent Belle data in the la level. The blue 
points are consistent with the recent Belle data on B — > rv in the la level [15]. The SM 
point of QV^sm = c r/ c sm = is slightly deviated from the red region, but it is in good 
agreement with the recent Belle data. If the new Belle data is combined together with the 
old data, the small discrepancy in BR(£> — > rv) might disappear. Then, BK(B — > tu) will 
give a strong constraint on the parameters related with charged Higgs boson. In Fig. 1 (b), 
we depict the allowed regions in our 2HDM for rrih+ and with both being scaled by 

tan j3. The red and blue regions are consistent with the combined data and new Belle data, 
respectively. 

In our U(l)' models, there exists an additional gauge boson Z'. One of the Higgs fields 
is charged under the extra ^(1)', so that the p parameter deviates from the SM prediction 
at tree level by 

A PtICC = {h^H^/v) 2 } 29 — 2 z 2 , (16) 

where m| = g 2 z v 2 and m|, = g' 2 v 2 {hf ((Hi) a/2 / v) 2 } + g' 2 h 2 v 2 . This form can be applied to 
2HDM, 3HDM and etc., fixing the charge assignment, {hi}, of Hi and of $. In the case 
of our 2HDM, hi((Hi)\/2/v) 2 = sin 2 /3. Therefore, the tree-level p parameter favors small 
tan/3 region [17]. This in turn implies that Y™ of 0(1) can be realized for mh+f tan/3 of 
O(100) GeV. 

On the other hand, the BaBar data on R(D) and R(D*) [4] prefers a large C c £ /C c g M (see 
Eqs. (14) and (15)). In Fig. 2 (a) and (b), we show favored regions (a) for Yukawa couplings 
\Y™\ and \Y™\, and (b) tan/3 and m h +, which are consistent with R(D) and R(D*) at 
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FIG. 1. Bounds on Yukawa couplings and m h + / tan (5 from B — > tv. The red region is the allowed 
region for the combined data [13] without the recent Belle data and the blue one is allowed by the 
recent Belle experimental result [15]. We used the relation YZ™ = V^V^Y^, ignoring the other 
elements of the Yukawa coupling. 



BaBar within la, respectively. The red points are consistent with the combined data for 
BR(Z? — > tv) while the blue points agree with the recent Belle data for BR(5 — > tv). \Y™\ 
is restricted to be less than 0.05 while \Y™\ is allowed to be 0(1). For the new Belle data 
\Y™\ is more constrained because the data is more consistent with the SM prediction and 
leaves little room for the charged Higgs contributions to B — > tv. In order to account for 
the discrepancies in R(D^), the Yukawa coupling \Y t ® u \ has to be sizable and its lower 
bound is about 0.2. As we have discussed in Fig. 1, \Y^\ of 0(1) might be consistent 
with BR(i? — > tv) experiments if R(D^) are not taken into account. Basically \Y™\ is 
constrained by BR(Z? — > tv) while \Y£ U \ is by R(D^), but they are related to each other 
through tan/3 and (g^) (see Eq. (5)). Hence, parameters to generate large Y™ are excluded 
by R(D^) data at BaBar, and the large top FB asymmetry at the Tevatron cannot be 
realized in our 2HDM. In Fig. 2 (b), tan > 3 is required in both the combined and new 
Belle data. For large tan/3, tends to be large. This is natural since C C ^ R in (14) and 
(15) are proportional to tan/3/m/ i + except for the last term in Eq. (14). 

Since only the RH up-type quarks are charged non-universally under t/(l)', our models 
do not generate the dangerous tree-level FCNCs contributions to B^-B^ B s -B s , and K°-K° 
mixings. The (u, c) elements of Yukawa couplings for neutral scalars and pseudoscalar, which 
may generate a tree-level FCNC contributing to D°-D° mixing, are small due to the sup- 
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20 40 60 80 100 



FIG. 2. Bounds on (a) \Y^ 1 \ and and (b) tan/3 and m h + in 2HDM. The points are consistent 

with R(D^) within la. The red points are consistent with the combined data of BR(Z3 — > tv) [13] 
while the blue points are in agreement with the new Belle data of BR(I? — > tv) [15]. 

pression factor of the light-quark mass. If a pseudoscalar has large (t, u) and (t, c) elements 
of Yukawa couplings, they may enhance a D°-D° mixing at the one-loop level. The loop 
of the pseudoscalar would induce the operators, C\(v^^cr)(^ur^cr) and C2(urCl)(urCl), 
but the contribution to C\ vanishes if external momenta are set to be zero. Therefore, only 
C*2 is non- vanishing, which is however suppressed by the factor, m 2 c /m 2 a . The upper bound 
on C 2 is discussed in Refs. [23, 24]: |C 2 | < 1.6 x 1CT 7 TeV -2 . Based on the general Higgs 
potential analysis in Appendix A, the mass difference between the charged Higgs and the 
pseudoscalar boson is at most the week boson mass scale. Roughly speaking, C<i could be 
estimated as 0((ml/m^)(Y t ( ^*Y t ^ u ) 2 / (16tt 2 )), so that we can expect that the points in Fig. 2 
do not disturb the SM prediction in D°-D° mixing very much. 

As discussed in Refs. [18], Y^ u ~ 0(1) is required to generate large Ap B and evade the 
strong constraint from the same-sign top-quark pair production at the LHC However, Fig. 2 
tells that Y™ should be less than 5 x O(10~ 2 ). Also, tan > 3 is required in the region 
with m h + > 200 GeV *, so that we cannot consider a sizable Z' interaction to achieve small 
Ap. Therefore in our 2HDM, it is difficult to find a favored region which is consistent 
with R(DW) and BR(B — > tv) as well as A^ B at the Tevatron. If R{D) and R(D*) become 
consistent with the SM prediction in the future, a tiny Y£ u is favored and the only constraint 
on Y^ will come from Fig. 1 (b). According to Fig. 2, small m h + / tan/3 is required by the 

•f The region m h + < 200 GeV is strongly constrained by B — > X s ^ and the exotic top decay. 
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large deviations of R(D^) and large Y£ u , so that a large / tan (3, where the new physics 
contribution is more suppressed, could be chosen to realize a large Y™, if the deviations in 
R(D*) become smaller. 



IV. 3HDM 

In order to enhance the Ap B , 0(1) (t, u) elements of Yukawa couplings for both a CP-even 
scalar and a pseudoscalar are required [18]. However, very tiny is required to explain 
the R(D^) in the BaBar experiments in our 2HDM. One simple solution to achieve both 
of the large Ap B and the BaBar discrepancies is to consider 3HDM. In the 3HDM of chiral 
U(l)' models, we have one more pair of charged Higgs and one more pseudoscalar, so that 
we can realize a scenario that one charged Higgs pair is constrained by B decays, but the 
other decouples with the B physics. In fact, we will find the parameter set that the Yukawa 
coupling of one pseudoscalar is 0(1). 



A. Yukawa couplings 

In this section, we consider the 3HDM with the (uk) = (1, 0, —1) U(l)' charge assignment 
of the RH up-type quarks, which was introduced in Refs. [18]. In our 3HDM, there are 4 
CP-even neutral scalars (hi, h%, h 3 , /i 4 ), 2 CP-odd pseudoscalars (ai, a 2 ), and 2 charged Higgs 
pair (hf, hf), after the gauge symmetry breaking. The Goldstone modes which are eaten by 
W and Z bosons are linear to the VEVs, ((Hi), (H 2 ), (H 3 )), and the orthogonal directions 
correspond to the mass eigenstates of pseudoscalars (£" 2 ) an d charged Higgs (£12)- Defining 
((Hi), (H 2 ), (H 3 )) = -j= (sin j3 cos 7, cos j3, sin j3 sin 7), we find that (£12) are given by 



^— cos cos 7^ 



4V = cosa a ,, 



V 



sin j3 
cos /3 sin 7 



+ sin a a . 



^ ' — cos f3 cos 7^ 



/-a.c 

4 2 = — sin a a 



sin /3 
cos f3 sin 7 



+ cos a a 



( sin 7 * 


y— cos 7 J 

sin 7 




^— cos 7 J 



(17) 



(18) 
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The mixing parameters, a ajC , relate to the terms, and Ay, in the Higgs potential, as 
described in Appendix A. The Yukawa couplings of pseudoscalars are given by 

Yl f X) = -Y^ cos a a + Yf } sin a a , (19) 

Y au(2) = yu(2) ^ + f u(l) (2Q) 

where and are defined as 

f « <iu v{^-i^w}' (21) 



The down-type quark sector and lepton sector Yukawa couplings are 



Y^ {1) = 5^ tan /3 cos cn , (23) 
// 

K""- 1 ^ = Sa — - tan B cos a a , (25) 



13 « v 

Y ad(2) = _ s rn^ ^ ^ 

J y 



Y$ m = -c%^tan/3sina a . (26) 



The magnitude of off-diagonal elements of is the same as the Yukawa coupling of 
2HDM, after replacing (Ri2R*2) by {gp>)ij- The Yukawa couplings of the charged Higgs are 
obtained from the relations (8) with replacement of a a by a c in 



B. R(D^) and BK(B tu) decay in 3HDM 

Now let us consider R(D^) and BR(i? — > tu) in our 3HDM. There are two charged 
Higgs pairs that contribute to C| 6 and C^ 6 , which are estimated at the charged Higgs mass 
scale as follows: 

^ - vm T t^B (Y*) [ ^ + ^1 + Y^^ 



1 '"2 ' \ "1 "2 



(27) 



cos a r sin a 



m b m T tan /3 ^ 1 ^ — . (28) 



°SM \ ""/).+ ""h\ 

R(D^) and BR(5 — > tv) depend on several parameters. In order to find the regions 
favored by BaBar and Belle data, we vary each parameter in the following range: 1 < 
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FIG. 3. (a) \Y t l u{2) \ vs. \Y t a u u{2) \ and (b) tan/3 vs. m h + for sina c = 1 in 3HDM. The red (blue) 
points in (a) are allowed points for 200 (400) GeV < m h + < 400 (1000) GeV. 

tan < 100, 200 GeV < m h + < 1 TeV, 200 GeV < m h + < 400 GeV, < a a , c < 2tt, and 
< \Ri2R%\, \Ri3R*is\ < 1, respectively, in addition to constraints on the Yukawa couplings, 



IV 



u(l,2) 



tu(tc) 



< 1.5. In the numerical analyses, we take tan 7 = 1 in order to realize a small tree- 



level contribution to the p-parameter §. After imposing the experimental constraints from 
R(D^) at BaBar [1], BR(B -> rv) [13], and the D°-W mixing, we could find parameter 
regions consistent with all the experimental constraints. For the bound on BK(B — > rv), 
we use the HFAG value, because the result does not change much even if we adopt the new 
data at Belle. 

Let us discuss a few specific cases for simplicity. First, we consider the case where 
cosa c = 1. In this case, C q ^ R /C q g M do not depend on m h +, so that decouples from 
B — >■ D^*'tv and B — > tv, but gives constraints on the parameters through the D°-D° 
mixing. Then, we can apply the same discussion as in the 2HDM by replacing (R i2 Rj 2 ) 
with (g^jij. In this case, |^ (1) (i? t2J R* 2 ) | = \Y™(R t2 R* u2 )\ is satisfied, so that |V^ (1) | should 
be smaller than 0.05 according to the discussion in our 2HDM. On the other hand, Y^ 2 ^ 
depends not only on (_R i2 i?* 3 ) but also on (i? t3 _R* 3 ) (see Eq. (22)). Since the combination 
{RizR*jz) certainly enhances Y^ , the heavier pseudo scalar boson a 2 exchanging could 
realize the large Ap B when cos a a is also 1 " . 

Secondly, we consider the case where sina c = 1. In this case, C^ L /C^ U is independent 

of hf. Then hf decouples from B — > D^ru and B — > tis, but gives constraints on the 

§ Note that Eq. (16) and {^((^) V2/v) 2 } oc (1 - tan 2 7) in the 3HDM. 
™ In the appendix A, the condition for cos a c = cos a a = 1 is discussed. 
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points in (a) are allowed points for 200 (400) GeV < m h + < 400 (1000) GeV. 

parameters through the D°-D° mixing. In Fig. 3, we depict |Y^ U( ' 2 ' ) | vs. |y^ 2 ' 1 | and tan/3 
vs. m h + for sina c = 1. In Fig. 3 (a), the red points are allowed points for 200 GeV < m h + < 
400 GeV, while the blue ones are for 400 GeV < m h + < 1 TeV, respectively. As we see in 
Fig. 3 (a), we find that small |V f "" ( ' 2 ^| as well as |V t " u( ' 2 ^| of 0(1) is in agreement with R(D^) 
and BR(Z? — > tv). As we have already discussed, (V^ 2 ^ of 0(1) might realize the large 
v4p B at the Tevatron. In this scenario, h\ contributes to the B decays, and the a2 of the pair 
of h\ could also enhance the Ap B . As we see in Fig. 3 (b), tan/3 tends to increase slightly as 
m h + becomes large, so that small tan /3 is also favored to increase the enhancement, because 
the mass difference between and a<i would be an order of electroweak scale. 

Finally, we consider the degenerate charged Higgs boson case with m h + = fn h +. In this 
case, Cf /Cf M and C q ^/Cf M are independent of a c and Y^ 2 \ We show the scattered plots 
for |V < " u( ' 2 ^| vs. |Y^"^| in Fig. 4 (a) and for tan (3 vs. m h + in Fig. 4 (b), respectively. 
We find that the tendency between tan (3 and m h + looks similar to that in 2HDM or in the 
cosa c = 1 case in 3HDM. As discussed in Refs. [18, 19], the small mass of the pseudoscalar, 
such as around 300GeV, is favored to enhance Ap B . As we see in Fig. 4 (a), one can find 
a few red points with 0(1) |V t " u< - 2 ^| and light m h + which could realize the large Ap B at the 
Tevatron. If m h + is heavier than 400 GeV, very large Xij in the Higgs potential are required 
as discussed in the appendix A. 

In conclusion, we can find the parameters regions which could explain the deviation of 
R(D^) at the BaBar experiment and be consistent with BR(I? — > tv) at the B factories 
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in 3HDM. We further investigates three interesting cases, (i) cosa c = 1, (ii) sina c = 1, and 
(iii) m h + = In all the three cases, we could find the allowed regions with large |, 
which in turn could generate a large top forward-backward asymmetry at the Tevatron. 



V. THE OTHER CONSTRAINT ON m h + 

It is well known that B — > X s ^y can give a stringent bound on the charged Higgs mass 
depending on the details of models. In type-II 2HDM, the lower bound is about 300 GeV 
< mh+ at next-to- next-to-leading order [25-28]. In our multi-Higgs doublet models with 
flavor-dependent U(l)' gauge interactions, the Yukawa couplings of the charged Higgs have 
extra parameters. The B — > X s ^ decay occurs through the loop diagram involving the top 
quark and the charged Higgs boson, where the relevant element of Yukawa couplings is the 
(b, t) element. According to Eq. (8), we can expect that the {b,t) element of the charged 
Higgs is governed by the (t, t) elements of the pseudoscalar bosons. In principle, the (t, t) 
elements of the pseudoscalar bosons have other mixing parameters, such as (g^)tt in 2HDM, 
which however are not directly constrained by the semileptonic and leptonic B decays. There 
is a theoretical relation, Kgfytql 2 = (.9R)qq(.9R)tt(<l = u,c), in our 2HDM, so that 0(1) (gfytc 
for #(£>(*>) requires 0(1) {g u R ) a . 

The bound on B — > X s ^ at leading order (LO) up to O((100 GeV / m h +) 2 ) is given by 

100GeV\\_ o , __, 2 1 /lOOGeVV ^ 

J V fn h + ) ~ 

(29) 

where two relations Y~ u = \/2V tb Y£ u and Y+ d = V2V tb Y^ b d = m b tan/3/w are used [29]. If 
we assume tan /3 = 1 and rrih+ = 300 GeV, then we obtain a constraint —0.077 < Y£ u < 
0.262. Therefore we can expect that (g%)tt can be 0(1) without conflict with the B — > X s ^ 
constraint. 



0.20 <{-[ 46.26 + 46.83 In ( — — J J Y£ u tan/3 + 9.00(Y t a t u ) 2 j I— J < 0.79 



VI. SUMMARY 



In this paper, we investigated the constraints from the semileptonic and leptonic B decays 
on our 2HDM and 3HDM, which were proposed in Refs. [18, 19] in order to accommodate 
the top forward-backward asymmetry (^4fb) observed at the Tevatron. In Refs. [18, 19], 
the U(l)'- charged extra Higgs doublets were introduced in order to generate the realistic 
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Yukawa couplings for the up-type quarks (especially the top quark mass) [18]. In the 
previous study, not only CP-even scalar bosons but also pseudoscalar bosons are required to 
have 0(1) (t, u) -element Yukawa couplings in order to achieve the large Ap B and to evade 
the strong bound from the same-sign top-quark pair production signal at the LHC [18, 19]. 
On the other hand, such a large (t, u) element of the pseudoscalar Yukawa coupling allows a 
large (b, u) element to appear in the charged Higgs Yukawa couplings. This implies that our 
model might predict large deviations from the SM predictions in B physics. For example, 
the (t, u) element of the pseudoscalar boson is indirectly constrained by the B — >■ tv decay. 

Recently, the BaBar collaboration reported the interesting results for R(D^) in the 
semileptonic B decays, B — > D^tv. The combined results deviate from the SM predictions 
by 3.4o" and require large FCNCs beyond MFV if the signal comes from charged Higgs 
exchanges. Our U(l) '-flavored multi Higgs doublet models naturally realize a large (b,c) 
element of the charged Higgs. Therefore in this paper, we investigated if our models explain 
both the BaBar discrepancies and A l FB at the Tevatron. 

In our 2HDM, only one Higgs doublet is charged under U(l)', so that small tan/3 is 
required to avoid the constraint on the tree- level p parameter. The R(D^) discrepancies 
at BaBar require large new physics effects, so that small m h + / tan /3 is required as we see 
in Fig. 2 (b). However, a small m h + / tan/3 requires a very small tan/3 according to 

the bound from B — > tv as shown in Fig. 1 (b), so that we could not find the points in 
Fig. 2 (a) with large which is needed for large Ap B . If the R(D^) converge to the SM 
prediction in the future, we would not need consider such a large new physics effect and a 
large Y£ u , so that we would be able to choose the points with a large and small tan/3 
in Fig. 1 (b). 

In our 3HDM, we have more freedoms: one more charged Higgs pair and one more 
pseudoscalar boson. It is not difficult to find the allowed points for R(D^) and BR(I? — > tv) 
in the general case. In the limit, cosa c — > 1, one of charged Higgs does not contribute to 
the (semi)leptonic B decay, so that we can describe the scenario that one of charged Higgs 
explains the BaBar results and the other becomes independent of the B physics. The 
explicit relation like Eq. (8) is not respected in this 3HDM because of the extra freedom, 
ce a ,c, but we also investigated the scalar potential in Appendix A and derived the condition 
for cosa a)C = 1, where both of Ap B at the Tevatron and R(D^) at BaBar could be achieved. 
In Figs. 3 and 4, we also plotted the allowed region, setting tan 7 = 1 which corresponds 
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to Aptrcc = 0, and discussed the two cases: sina c = 1 and m h + = m h +. We can find the 
points with 0(1) (t,u) Yukawa coupling and 300 — 400 GeV charged Higgs mass, which 
are consistent with R(D^) and D°-D°. tan/3 < 40 should be satisfied in sina c = 1 case, 
and the pair, (^2,02), could achieve the BaBar discrepancy and Ap B . In the degenerate 
limit, m h + = m h +, the a c dependence disappears in the B decays and we find the tendency 
between tan (3 and m h + as shown in Fig. 4. However, the new contributions from the extra 
scalars are also constrained by the D°-D° mixing, so that only a few points with a large 
(t, u) Yukawa coupling are allowed if m h + is lighter than 400 GeV. 

We also commented on the other bound on the charged Higgs boson. One of the most 
important bounds on the charged Higgs boson comes from B — > X s ~f. In our models, the 
(t, t) element of Yukawa coupling contributes to the process at LO and we derived the bound 
in the limit that the charged Higgs boson is heavy. As discussed in Refs. [25-28], the lower 
bound on the charged Higgs mass in the type-II 2HDM is around 300 GeV, but the mass 
below the bound is possible in our models, depending on the {t,t) element. 

Before closing this paper, let us emphasize once more the importance to study phe- 
nomenology in the framework of a well defined consistent renormalizable lagrangian. The 
original Z' model for the top forward-backward asymmetry has been excluded several times 
by the upper bound on the same-sign top-pair production cross section. However the model 
with only Z' is not well defined since it is not renormalizable, and not realistic because the 
up-type quarks including top quark are massless. It is mandatory to extend the Higgs sector, 
by introducing new Higgs doublets with nonzero U(l)' charges. Making such an extension 
actually affects the top phenomenology a lot. Basically all the top-related observables are 
affected by extra Higgs doublets, both neutral and charged Higgs bosons. Also the B meson 
and D meson sectors are modified too. Maybe it is timely to remind ourselves that the SM 
with three quarks u, d, s produces too large FCNCs in the kaon sector, and is immediately 
excluded. This disaster can be cured only by introducing the 4th quark, the charm quark, 
which also solves the problem of gauge anomaly. The experience with the charm quark 
already teaches us that it is important to work in a minimal consistent (anomaly free) and 
renormalizable model in order to do a meaning phenomenology. 
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Appendix A: Potential Analysis 

The general potential in multi-Higgs-doublet models with U(l)' is given by 

V = m\R\E { - [m%^)E\U 3 + h.c.) 

+^(HjH, i ) 2 + Xij(HjHi)(HjHj) + Xij (H}Hj) , (Al) 

where m^($) = = Xa = 0, A^- = Xji and A^- = X^ are satisfied. m?-($) is the function 
of $, which is the complex scalar to break U(l)' and the function is fixed by the charge 
assignment of the fields. Each Hi includes neutral, pseudoscalars and charged Higgs, 



.72 ^ 72* 



Vi satisfies the stationary condition, 



vf , V,jV 2 ~ ViV 2 



= mfvi - (ml + m 2 *) Vj + A^ + + (A3) 

Assuming $ gets nonzero VEV, the mass matrices of pseudoscalars and charged Higgs are 

(M a % = -Mj. + ^M^, (A4) 

(M 2 h+ h = (M a % + X^-f- - \ ik V l5 ih (A5) 

where M% = m% + rnBl and we assume that Mf 4 is real. 
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In 2HDM, we can find one simple relation between m 2 H and m 2 , 



m 



H+ 



^12 — - 



(A6) 



This means that the mass difference is at most the electroweak scale. 

In 3HDM, the directions of massive modes in pseudo and charged Higgs sectors can not 
be fixed, as we discuss in Sec. IV. The condition for cosa c>a = 1 is obtained from the 
following calculation, 



sin 7 » 



Mi 



\ 







cos 



V 



sin 7 * 



Ml, 



\ 





cos 7 J 



( 



t^tan,o 1 Jcos7 

— M 2 2 sin 7 + M| 3 cos 7 
M? 1 



\ 



M 



-^1 2 3 — 

J-O SII17 



sin 7 » 





cos 7 J 



X 



23 tan /3 tan 7 
/ 

2 



V 

y 



sin 7(Ai2 cos 2 + Ai 3 sin 2 (3) 



\ 



(A 2 3 — Ai 2 ) cos (3 sin (3 cos 7 sin 7 
^ - cos 7(A 2 3 cos 2 (3 + A 13 sin 2 (3) j 



(A7) 



(A8) 



That is, M 2 2 tan7 = M| 3 cos 7 and A23 = A12 are required. 
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